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Lewis number and Damkohler number 
effects in vortex-flame interactions 

By J.-M. Samaniego 
1. Motivation and objectives 

Premixed flames are encountered in numerous practical configurations, and a 
better knowledge of the fundamental mechanisms controlling their propagation is 
needed to help model their behavior. In practical devices, premixed flames gener- 
ally propagate into a turbulent stream of reactants, and their structure and local 
flame speed are modified. Although some recent studies demonstrate the persis- 
tence of laminar flame structures under highly turbulent conditions (Poinsot et ai 
1991, Furukawa et ai 1993), the structure of turbulent premixed flames, either 
as a flamelet or thickened flame, is still a matter of debate. Hence, the develop- 
ment of turbulent combustion diagrams is an ongoing effort (Borghi 1988, Poinsot 
et ai 1991). However, despite uncertainties concerning the actual structure of pre- 
mixed flames in real combustion devices, studies of the propagation properties of 
a laminar flame in turbulent flow conditions contribute to a better understanding 
of turbulent premixed flames. In such flows, premixed flames experience stretch, 
strain, curvature, and unsteadiness, and are subject to various effects depending on 
the thermodiffusive properties of the mixture (Lewis number), adiabaticity (heat 
losses), and detailed kinetics (Clavin 1985, Borghi 1988, Law 1988, Clavin 1994). 

In order to study these effects in detail, a combined experimental-numerical study 
of the interaction of a two-dimensional vortex pair with a plane premixed laminar 
flame has been carried out. The selected geometry is two-dimensional in order 
to allow for the use of quantitative line-of-sight measurement techniques and for 
comparisons with two-dimensional direct numerical simulations. The experiment is 
used to identify possible effects of the Lewis number and of radiative heat losses, and 
direct numerical simulations reproducing the experimental conditions are used to 
investigate the role of the Lewis number, of heat losses, and of multi-step kinetics. 
The heat losses considered in this study are essentially radiative losses from the 
burnt gases and do not include conductive losses to the walls. 

This report presents the final results of this study of vortex-flame interactions. 
The main result is the creation of experimental data sets of vortex-flame interac- 
tions, which demonstrate the importance of the Lewis and Damkohler numbers to 
the flame and which allow for comparison with direct numerical simulations. In turn, 
the comparison between experimental and numerical results, presented in Mantel 
1994, shows that the Lewis number effects and complex chemistry play significant 
roles, which can be reproduced using a two-step reaction mechanism. 
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FIGURE 1. Schematic view of the facility 

2. Accomplishments 

2.1 Experimental details 

An experimental facility with a two-dimensional flow has been developed and has 
been described in the last issue of the Annual Research Briefs, where details about 
the vortex generation and flame stabilization are reported (Samaniego 1993). The 
two-dimensionality of the flow field is Em important feature of this study since it 
allows for comparison with two-dimensional direct numerical simulations of vortex- 
flEune interactions. The choice of a two-dimensional flow field is dictated by the 
necessity of providing numerous quantitative experimental and numerical data in 
the simplest though relevant flow configuration by the use of line-of-sight COJ 
emission imaging eis a qusmtitative measurement of the heat release rate, and by 
limitations of current computer capabilities that would render three-dimensional 
calculations with multi-step kinetics too expensive. 

The test section comprises a vertical duct, with a square cross-section of 63.5 x 
63.5 mm, equipped with quartz windows for optical access (see Fig. 1). Mixtures 
of fuel and oxidizer are fed into the test section through a contoured converging 
nozzle. Different fuel and oxidizer mixtures Eire used to investigate the effect of the 
Lewis number. The Lewis number is determined using mixture-averaged properties 
(Strehlow 1979, Williams 1985). Five fuel-lean mixtures were investigated: <j> — 0.55 
C H a / Air (Le = 0.96), <f> = 0.46 C 2 H 4 /Air (Le = 1.24), <f> = 0.51 C 3 H S /Air (Le = 
1.60), <f> = 0.61 CH 4 /0 2 /C0 2 /Ar ( Le = 0.80), and <j> = 0.61 CH 4 /0 2 IC0 2 j He 
(Le = 1.54). The flame is stabilized on a heated Nichrome wire of 0.5 mm diameter, 
resulting in a V-shaped flame. 

Several parameters describing the vortex pair suid the flame control the interaction 
between the vortex and the flame. The vortex parameters, determined from smoke 
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Table 1. Vortex parameters 


type of vortex 

s (cm) 

Vo (cm/s) 

R (cm) 

r r (ms) 

Slow 

0.85 

135. 

10.0 

4.15 

Medium 

0.60 

350. 

4.0 

2.07 

Fast 

0.60 

600. 

7.4 

1.09 


Table 2. Flame parameters 


Mixture 

<t> 

Tad 

(A) 

(deg) ( 

Sl 

cm/s) 

Sl 

(cm) 

Le 

HL 

xlO -3 

C Hi/ Air 

0.55 

1567 

13 

9.5 

0.022 

0.94 

0.55 

C 2 H 4 /Air 

0.46 

1515 


9 

0.022 

1.25 


C 3 H S /Air 

0.51 

1530 

15 

9 

0.020 

1.60 

0.40 

CHi/0 2 /C0 2 /He 

0.61 

1643 


13.8 

0.033 

1.54 

1.10 

CH 4 /0 2 /C0 2 /Ar 

0.61 

1615 


6.7 

0.024 

0.82 

2.50 


visualization experiments, are: the distance between the vortex centers, s, the vortex 
pair self-induced velocity, Vd> and the radius of curvature of the trajectory, R. The 
flame parameters are: the mean flame angle, a m , the flame speed, Sj,, the flame 
thickness, 6l, the Lewis number, Xe, and a heat loss coefficient, HL , which is 
defined as the ratio of the radiative loss in the reaction zone to the chemical energy 
release. They depend on the mixture composition (fuel, diluent) and on the bulk 
flow velocity, Vo- In the present study, Vo = 0.35 m/s. The vortex and flame 
parameters are reported in Tables 1 and 2. Details on the determination of these 
parameters can be found in Samaniego et al 1994a. 

Line-of-sight CO\ emission imaging, using an intensified CCD camera, is per- 
formed to obtain heat release rate fields. Each image is digitized in a 184 x 240 
8-bits pixel array (a pixel corresponds to a field of vue of 0.364 x 0.364 mm). The 
light intensifier is used as a fast shutter in order to freeze the flow field. The expo- 
sure time is 250 /is, resulting in a spatial resolution of 1.5 mm in the worst case. A 
timing circuit allows the triggering of the intensifier at selected instants during the 
interaction. For each value of the time delay, ten to fifty realizations are averaged 
in order to improve the signal-to-noise ratio. Images are corrected for background 
noise and pixel-to-pixel variations. After correction, the uncertainty on the pixel 
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value is ± 2 % of the maximum value in the worst case. 


2.2 Determination of the heat release rate 

The rate of heat release is an important flame quantity which controls its prop- 
agation, and its monitoring provides insight into the respective effects of the Lewis 
number, radiative losses, and complex chemistry. The choice for studying this 
quantity is made for two reasons: first, it can be determined experimentally by 
measuring the flame chemiluminescence; second, this quantity is predicted by most 
combustion models and can therefore be used for their validation by comparison 
with experimental data. 

Chemiluminescence refers to light spontaneously emitted by electronically excited 
molecules as they decay to their ground state. In hydrocarbon flames such as those 
investigated here, flame chemiluminescence occurs in the UV and visible range of 
the spectrum and is due principally to four emitters: CH*, C OH* , and C0\. 
These electronically-excited molecules have different emission spectra. The first 
three emitters have band- structured emission spectra which correspond to specific 
electronic transitions ( CH * peak at 431.5 nm corresponding to the A 2 A — ► X 2 IT 
transition, C% peak at 516.6 nm corresponding to the A 3 U — X 3 Tl transition, OH* 
peak at 306.4 nm corresponding to the 2 E + — > 2 II transition Gaydon 1974), while 
COX has a broadband emission spectrum extending from 250 to 800 nm (Mvers &: 
Bartle 1967). 

Recently, Samaniego et al. 1994b, after showing that CO\ is the main emitter 
in fuel-lean flames, have successfully correlated the COJ emission intensity and the 
rate of heat release for a <j> — 0.55 CH\/Air flame based on complex chemistry 
calculations of strained flames. The resulting correlation takes the form: 



where Q is the integral of the heat release rate, across the flame front, i.e.: 



I is the CO\ emission intensity integrated across the flame front, i.e.: 


1 = 



and the subscript 0 refers to the unstrained premixed laminar CH^/Air flame (<f> = 
0.55). 

It was shown that expression ( 1 ) is not affected by unsteadiness and can be used to 
determine the heat release rate in situations where strain-rate is the only significant 
parameter controlling the flame evolution. As is discussed later, it is the case in the 
present study of vortex-flame interactions since the effect of curvature, which is not 
accounted for in Eq. (1), can be neglected. 



Experiment on vortex-flame interactions 


33 





FIGURE 2. Dithered images of CO $ emission of a <f> = 0.51 propane-air flame, left: 
t = 0, center: t = 17 ms, right : t — 19 ms. 

However, Eq. (1) is only strictly valid for the <j> = 0.55 CH^/Air flame, and new 
correlations should be sought for other mixture compositions (i.e. different fuel, 
different diluent). 

2.3 CO\ emission imaging 

Instantaneous images of CO £ emission from the flame were obtained using an 
intensified Amperex CCD camera equipped with a NIKON glass lens (/ = 105 
mm, //1.2) and a band-pass filter BG39 with cut-off wavelengths at 340 and 700 
nm. This filter is used to reject near-infrared emission (700-1200 nm) from the 
burnt gases originating from H^O excited molecules (Gay don 1974). 

Fig. 2 shows a sequence of three flame emission images of a C^Hs/Air flame at <f> 
= 0.51 taken during a vortex-flame interaction. The parameters of this interaction 
are: s — 6.0 mm, Vb = 3.5 m/s, R = 40 mm, Si ~ 0.10 m/s, 6i ~ 1.0 mm, Le 
= 1.60, a ~ 15 degrees. 

The first image shows the unperturbed flame fronts. The next two images in 
Fig. 2 are taken respectively at t = 17 ms and t = 19 ms. Three regions in the 
distorted portion of the flame can be identified: 1) a region of extensive strain-rate 
and negative curvature, located just ahead of the vortex pair. In this region, strain- 
rate effects are predominant, and Eq. (1) can be confidently used to infer the heat 
release rate; 2) &: 3) regions of high positive curvature, located above and under 
region 1. In these regions, curvature may play a significant role, and the validity 
of Eq. (1) remains to be evaluated. Therefore, quantitative comparisons between 
experimental and numerical data should be restricted to the regions of high strain. 

For this purpose, X mgn is defined as the minimum value of Z(<r), where Z(<r) is 
defined as the integral of the CO\ emission, /, across the reaction zone, at the 
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FIGURE 3. Slow vortex. Evolution of E. o : Le = 0.94, & : Le = 1.25, • : 
Le = 1.60 (see Table 2). 

location a along the flame front, i.e.: 



where J_ indicates that integration is performed along a normal-to-the-flame front, 
and 


Imin = min 1(a) 

O 

In all cases, X mi „ occurs in front of the vortex pair in the region of high extensive 
strain. The corresponding minimum heat release rate, Q m i n is obtained from X mtn 
using expression (1). 

2.4 Effect of the Lewis number 

The theory of stretched flames predicts that the behavior of strained flames is con- 
trolled by the Lewis number, Le (Clavin 1985, Law 1988). For positively stretched 
flames such as in vortex-flame interactions, the burning rate should increase for 
Le < 1 flames while it should decrease for Le > 1 flames. This prediction is checked 
by monitoring the evolution of I m ; n , during vortex-flame interactions, between the 
slow vortex and all flames (each flame corresponds to a different Lewis number - 
see Tables 1 and 2). 

It is first verified that all flames undergo the same stretch history. Fig. 3 shows 
the evolution of the flame surface area, E, during vortex-flame interactions with the 
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FIGURE 4. Slow vortex. Evolution of K m ean • o : Le — 0.94, © : Le = 1.25, • : 
Le = 1.60 (see Table 2). 



0 5 10 15 20 25 30 

time (ms) 

FIGURE 5. Slow vortex. Evolution of l mxn for all investigated flames (see Table 
2). □ : Le = 0.82, o : Le = 0.94, « : Le = 1.25, e : Le = 1.54, • : Le = 1.60. 


slow vortex pair for all flames. The flame surface area, E, starts increasing after t 
= 15 ms. A mean stretch-rate, A' mean , can be derived: 
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l_dZ 
E dt 


i\ mean also starts increasing after t = 15 ms and reaches about 100 to 150 s~ x 
at t = 25 ms (see Fig. 4). All flames experience the same overall stretch history 
which allows us to interpret possible differences between these flames in terms of a 
Lewis number effect. 

The effect of the Lewis number on the evolution of J m i n for all flames is shown 
in Fig. 5. In all cases, I m i n starts at the value of the unperturbed laminar flame 
and then decreases under the action of the vortex pair. It can be noticed that the 
decrease is more pronounced for higher Lewis number flames. This is better demon- 
strated in Fig. 6 where the maximum slope of the evolution of I min , 6, is plotted as 
a function of the Lewis number, where the slope b decreases with increasing Lewis 
number. 



Le 


FIGURE 6. Slow vortex. Correlation between the rate of variation of I m ,„ (6), 
with the Lewis number (Le). 

These results are only in partial agreement with the asymptotic theory: higher 
Lewis numbers lead to flames that are weaker when strained. However, the theory 
predicts that the burning rate of Le < 1 flames increases with strain, but we observe 
a decrease of J m ; n for all flames, including the Le = 0.8 flame. This may be 
attributed to various causes: unsteady effects (see section 2.5, theory limited to low 
stretch-rates, complex chemistry effects (see DNS). 

2.5 Effect of the Damkohler number 

Unsteadiness is an important factor in these vortex-flame interactions which is not 
accounted for in the theory of stretched flames. Unsteadiness is quantified by the 
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Damkohler number, defined as the ratio of a mechanical time, r m , over a chemical 
time, r c , i.e. 


r c 

In this study, r m is the characteristic time of evolution of the flame surface area 
and r c is the residence time through the flame, i.e. r c = 61 / Sl- T m is defined by 
the time for E to increase from 1.25 x E 0 to 2.00 x E 0 where E 0 is the unperturbed 
value. 

2.5.1 Observations 

The effect of the Damkohler number is determined by investigating the evolution 
of the overall heat release rate Q toi , and the minimum heat release rate Q m ,„. Qtot 
is defined as the total heat released by the flame of surface area E and is computed 
as: 

Qtot = J Q{cr)da (2) 

Fig. 7 shows the evolution of E, Qtot, and Q m in for the <f> = 0.55 CH^/Air flame 
during its interaction with the slow, medium, and fast vortex pairs (see Tables 1 and 
2). These interactions are characterized by Da = 1.80, 0.90, and 0.47, respectively. 
Lower values of Da can be interpreted as higher stretch rates. Each set of curves 
corresponds to a different Damkohler number. In all cases Q to t lags E. This is 
consistent with the fact that Q m i n decreases due to flame stretch. Indeed, if Q(cr) 
were constant along the flame front during vortex-flame interactions, Qtot would be 
proportional to E (see expression (2)). Hence, the non-linearity of the relationship 
between Qt ot and E is a consequence of stretching due to the vortex pair. 

In order to stress the effect of Da, E and Qtot are plotted in a non-dimensional 
time- frame, where time is non-dimensionalized by r m , and the origin corresponds 

to the time at which E = 1.25 x E 0 (Fig. 8). The evolution of E is independent 
of the Damkohler number. This follows mainly from the non-dimensionalization 
procedure which imposes E(£ + = 0) = 1.25 x Eo and E(f*~ = 1) = 2.00 x Eo. 
The evolution of Q to t shows a certain similarity in all three cases: Qtot starts 
decreasing, reaches a minimum value, and finally increases at a pace comparable 
with the evolution of E. One can notice that the minimum value reached by Q to t 
decreases with decreasing Da and that higher values of Da lead to less non-linearity 
between Q to t and E. Third, the decrease Q m in is less pronounced for higher values 
of Da. These observations display trends that are in agreement with the idealized 
infinitely fast chemistry limit, or Da — ► 00 , where the internal structure of the flame 
is not affected by the flow. In this case, Q(a) would be constant, and Q toi would 
be proportional to E. Indeed, as shown in Fig. 8, Q to t converges towards E as Da 


increases. 



38 


J.-M. Samani 



Figure 7. CH^/Air flame ( <j> — 0.55). E\ 
(• ), and Qmin / Qo (• )• Top: slow vortex ( Di 
(Da = 0.90), bottom: fast vortex (Da — 0.47). 
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Figure 8. CH^/Air flame (4> = 0.55). Effect of Da on the relationship between 
E/E 0 and Qtot/Qtot, o- Qtot/Qtot, o -» o : Da = 1.80, « : Da = 0.90, • : Da = 0.47. 
E/E 0 — > □ : Da = 1.80, a : Da = 0.90, ■ : E>a = 0.47. 


2.5.2 Stretching-relaxation mechanism 

For finite values of Da, unsteadiness, in addition to stretch, plays an important 
role. This can be seen in recent studies of unsteady strained laminar flames (see for 
example Darabiha 1992, Egolfopoulos 1994). A simple stretching-relaxation model 
accounting for unsteadiness may be derived. The main assumption is that two 
mechanisms are important: unsteady strain which tends to thin the flame without 
modifying the internal flame structure, and the diffusion-reaction process which 
relaxes the flame back to an unstrained laminar structure. Effects of curvature are 
neglected. Both mechanisms contribute to dQ/dt. They can be modeled by two 
additive effects: a term of the form — KQ due to stretch, and a relaxation term of 
the form (Q 0 — Q)/r r , where r r is a relaxation time. The term due to strain-rate 
can be justified by two arguments: 1) curvature being neglected, stretch is equal 
to the strain-rate tangential to the flame front, and we have K = ax', 2) the flame 
front undergoes a compression —ax along its normal. Therefore the flame front 
is compressed at a rate -K along its normal, and the rate of change of Q due 
to stretch is equal to —KQ. The second term is an ad-hoc representation of the 
diffusion-reaction process using a linear relaxation model and is justified by the fact 
that the Q must relax to Q 0 in the absence of stretch (actually, different forms for 
the relaxation term could be used since there is no real justification for a linear 
relaxation process). Therefore, an evolution equation for Q may be derived: 


dQ 

dt 


-KQ + 


Qo — Q 


T i 


( 3 ) 
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FIGURE 9. Evolution of Q/Qo with strain-rate for a steady strained CH^/Air 
flame (<f> = 0.55). Comparison between model predictions (with r r = 1.25 mi) 

and complex chemistry calculations. : complex chemistry, : linear 

relaxation model, : quadratic relaxation model, : cubic relaxation 

model. 

Interestingly enough, in the steady case (i.e. ^ = 0), this model yields: 

Q = i 

Qo 1 + Kr r 

This prediction is compared to the evolution of the relative heat release rate for a 
strained <j> = 0.55 CH^/Air flame, computed from complex chemistry calculations 
(Samaniego et al. 1994b) (Fig. 9). Other predictions using non-linear relaxation 
terms are also plotted for comparison (quadratic: ^°(1 — Q/Qo) 1 ^ 2 , and cubic: 

^(1 — Q/Qo ) 1 ^ 3 )* In this case, r r = 1.25 ms is the same for all three relaxation 
models and is chosen to best reproduce the evolution of Q/Qo predicted by the 
complex chemistry calculations. It appears that all relaxation models predict cor- 
rectly the qualitative behavior of Q/Qo and that quadratic relaxation provides the 
best quantitative agreement. In the following, the linear relaxation model, although 
imperfect, is used to further investigate the relationship between Q\ 0 t and E. 

Eq. (3) may be integrated over E using a transport theorem which accounts for 
the fact that E also changes in time (see Candel & Poinsot 1990). One obtains, 
after proper non-dimensionalization: 

£*£+<«.-* ( 4 ) 

where t + = t/T m , C = r r /r c is a constant, Q~l ot = Q to t/Qtot 0 » and E + = S/E 0 . 
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Figure 10. Evolution of Qtot I Qtot, o during vortex-flame interactions. Compar- 
ison between the linear model (r r = 1.25 ms) and measurements. Model (Eqs. 5 

and 6) : S+, : Qtot for Da = 1.80, : Qf ot for Da = 0.90, 

: Qf ot for Da = 0.47. Experiment — ► e : E + , o : Qf ot for Da = 1.80, « : Qf ot 

for Da = 0.90, • : Qtot for Da = 0.47. 

From this expression it can be seen that, in the limit Da — > oo, Qtot — E"*" and the 
total heat release rate is proportional to the flame surface area, and, in the limit 
Da —* 0, Qt 0 t = 1 and the total heat release rate is constant. 

In turn, Eq. (4) can be integrated once the evolution of E is known. A convenient 
analytic expression for is provided by: 


E + = 1 + 4 t_1 (5) 

This formulation mimics the exponential-like growth of E + in Fig. 8, and satisfies 
E + (0) = 1.25 and E + (l) = 2. Integrating (4) yields: 


Qtot = i + 


Da 

(Da + C In 4)' 


it-i 


(6) 


Fig. 10 shows a prediction of the effect of the Damkohler number using these 
analytic expressions for E and Q to t, with C = 0.25. Despite some discrepancies, in 
particular, the inability of the model to reproduce the initial decrease of Qtot, the 
effect of Da on Q to t is qualitatively well reproduced. 


3. Conclusion 

This experimental study of vortex-flame interactions has been successfully con- 
ducted, and the initial goals of this study have been achieved. The results have 
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demonstrated the importance of the thermodiffusive diffusive properties (Lewis 
number) and of unsteadiness (Damkohler number) on the response of the flame. 
Quantitative comparisons with direct numerical simulations have demonstrated the 
importance of complex chemistry effects which could be reproduced using a two-step 
mechanism. 

Further experimental studies of vortex-flame interactions could be carried out to 
elucidate the effect of the flame (heat release) on the vorticity field. Preliminary 
results show that the vortex pair is significantly altered in the presence of a flame, 
and this might be due to the baroclinic torque as well as to buoyancy. This study 
would shed new light upon the interaction between turbulence and heat release, 
such as possible jump conditions for turbulence statistics through the flame. 
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Erratum 

The author would like to point out two errors that were made in the previous 
issue of the Annual Research Briefs (Samaniego 1993). First, what he mistakingly 
referred to as CH* emission is actually COJ emission. Second, the presumed linear 
relationship between the so- denominated CH * emission and reaction rate does not 
apply. Both these errors were made before the study on COJ emission, which 
identified COJ as the main emitter in these flames and led to a quantitative (non- 
linear) relationship between COJ emission and heat release rate (Samaniego et al. 
1994b). 
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